Coherence of a Josephson phase qubit under partial-collapse measurement 
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We discuss quantum evolution of a decaying state in relation to a recent experiment of Katz et 
al. Based on exact analytical and numerical solutions of a simple model, we identify a regime where 
qubit retains coherence over a finite time interval independently of the rates of three competing 
decoherence processes. In this regime, the quantum decay process can be continuously monitored 
via a "weak" measurement without affecting the qubit coherence. 
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An impressive recent progress in implementing sim- 
ple quantum systems related to quantum computing]!]] 
permits to address experimentally the long-standing con- 
troversies on quantum measurement^], j3|. In particu- 
lar, it is becoming possible to study directly what hap- 
pens "inside" the quantum state collapse during a con- 
tinuous (weak, partial, incomplete, etc.) measurement. 
While continuous measurement (CM) of an ensemble of 
quantum systems can be described just as an ensemble 
decoherence H Q, a much more subtle and interesting 
topic is CM of a single quantum system. 

This topic is well-developed in quantum optics j^], with 
the most advanced experiments including a demonstra- 
tion of a continuous quantum feedback ||. While the 
formalisms and terminology used by different groups in 
relation to the continuous (partial, etc.) collapse are 
quite diverse, the most widely known theoretical ap- 
proaches are so-called POVM (positive operator-valued 
measure^]) and "quantum trajectory" [||. In condensed 
matter physics a similar approach has been introduced 
as the "quantum Bayesian" formalism[Q. 

The first direct condensed matter experiment on par- 
tial collapse has been realized recently^]. (Somewhat 
similar experiment was proposed in optics) 3 but never 
realized.) The experimental setup of Ref. 8 was based 
on the Josephson phase qubit ]lCfl [Fig. 1(a)], which has 
an asymmetric double-well potential profile. Two lowest 
levels (with energies Eo and E\ ) in the shallow "left" well 
were used as qubit states |0) and |1) [Fig. |l|(b)]. The lev- 
els in the deep "right" well were significantly broadened, 
essentially creating a continuum of states. With some 
(over)simplification, the experiment can be presented in 
the following way. The qubit was prepared in a super- 
position state V(0) = cvo(0)|0) + ai(0)|l), and then the 
barrier was lowered for a time t to allow a partial tunnel- 
ing from the state |1) into the right well (Tt ~ 1, where 
F is the tunneling rate). Selecting only the cases when 
the tunneling had not happened ( "null-result" ) , the qubit 
state was subsequently examined by the quantum-state 
tomography. Experimental results S were consistent 111] 
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FIG. 1: (a) Schematic of the phase qubit: superconducting 
loop interrupted by a Josephson junction, (b) Potential profile 
and level structure of the phase qubit. Upper level |1) decays 
with the rate Y. (c) Idealized level structure used here. Left- 
well levels with energies Eo and E\ form a qubit. The state |1) 
can decay into the right well where the average level spacing 
is A and the bandwidth is A. The tunneling amplitude to the 
state \n) is f n . 

with the simple formula 

m = MO) e- iEot \0) + ai(0) e-^ f e- rt / 2 |l) 
v /|a (0)| 2 + |a 1 (0)| 2 e- r * 

which follows from the quantum Bayes rule@, ]l2) for par- 
tial measurement of the qubit. Notice that for Tt 3> 1 
this formula describes the "orthodox" projective collapse 
onto state |0) (this regime is usually used for the phase 
qubit measurement! 10 by sensing the tunneling into the 
right well with a nearby SQUID), while for Tt ~ 1 the 
collapse is only partial. Therefore, the experiment || has 
shown that after the partial collapse the qubit remains 
almost perfectly pure, while its evolution is information- 
related; in particular, the amplitude of state |0) gradually 
grows without "physical" interaction. 

The purpose of this work is to understand why and how 
well a metastable qubit may retain coherence despite de- 
cohering processes in its environment. This is important 
for understanding of the partial-collapse measurement in 
Ref. |^, but even more so for future experiments on con- 
tinuous monitoring for qubit decay where measurement- 
induced decoherence would be inherent. We focus on a 
simplified model with the level structure as in Fig. 0(c), 
where the (qubit) states in the left well experience no 
direct decoherence, whereas those in the right well and 
the tunneling Hamiltonian are subject to decoherence. 
As we discuss, in this case the qubit remains pure as 
long as the tunneling out of the left well is an irreversible 
process. We use analytical and numerical techniques to 
illustrate situations where such an irreversibility is due 
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to the choice of system parameters (e.g., for nearly con- 
tinuous spectrum in the right well) or where it happens 
dynamically due to the evolution properties in the right 
well. 

We write the system Hamiltonian in the block form, 



H 



H L T 
Tt H R 



(2) 



Here Hl is the two-level Hamiltonian in the left well, 
Hl = diag{i?o, -Ei}, Hr is the Hamiltonian in the right 
well, Hr = J2 n E n \n)(n\, and T is the corresponding tun- 
neling Hamiltonian (here and below £ denotes summa- 
tion over the right- well states only). Unless mentioned 
otherwise, we assume that only the transitions from the 
upper qubit state are allowed, T\ n = f n ^ 0, while the 
state |0) is fully disconnected, Tb„ = 0. 

Let us start with the simplest case of no decoherence 
during the tunneling time interval t, followed by an ideal 
orthodox quantum measurement which distinguishes left 
and right wells (technical realization of such measure- 
ment is discussed below). Then the state of the system 
can be described by a wavefunction -0 = (ipL, ^r) \ which 
starts as t/>i(0) = a (0)|0) + ai(0)|l) and ijj R (0) = 0, 
evolves according to the Schrodinger equation ip = —iHip 
before the measurement, and finally undergoes orthodox 
projective collapse at time t. In particular, the left- 
well component ipL is either zeroed if the escape is de- 
tected, or is rescaled to become the new wavefunction 
of the system]^, [l3| if the measurement finds no escape: 
(V'isV'ii) ~* (V'l/II^lIIi 0). In this case it is trivial to see 
that the qubit remains fully coherent in the interesting 
for us null-result scenario of no escape. 

Before the measurement the left-well components cyq, 
ct\ evolve as ao(t) = e~ %Eot ao(0) and 



£*i(t) 



iai(0) 
2tt 



de e 



e-^l -En\fn\ 2 /(t-E n +i0) : 
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In the intermediate range, A -1 <C t <C A -1 , the level dis- 
creteness is unimportant and the residue summation can 
be approximated by an integration. We obtain 



ftl (() = z ai (o)^ i+ffii V r '/ 2 , 



(4) 



where the decay rate is T = 2irD(e) 1 while the energy 
shift 8Ei and the prefactor Z are given by the integrals 



8E l =V 



dED(E) 
e-E ' 



1 

— = 1 

Z 



dED(E) 
(e-E + iO)* 



all evaluated at e = E\ + 8E\ — iT/2. Here we in- 
troduced the smoothened tunneling density of states 
(TDOS) D(E) instead of D(E) = \f n \ 2 S{E - E n ); 
unlike in Ref. |l4]we assume D(E) to have no discontinu- 
ities or singularities. 

Replacing D by / 2 /A, the decay rate can be written 
as T = 27r/ 2 /A, and the evolution is well-exponential 
only if A <C / <C V AA. In this case (which is essentially 
tunneling into continuum) we obtain the simple formula 
(|l|) with small corrections 8E\ and Z. 

Now let us add decoherence processes into the picture. 
We will consider only decoherence in the right well and 
between the wells, excluding explicit left- well decoher- 
ence which has a trivial effect. 

A simple model describing decoherence of right-well 
levels can be introduced by adding imaginary parts to 
their energies, E n — > E n — ij n /2. Physically, this cor- 
responds to processes of energy relaxation to additional 
levels which do not interact with |1). Then the wavefunc- 
tion formalism [Eq. (^)] is still valid, so the qubit remains 
pure after the null-result measurement, while the con- 
ditions for the exponential decay are now more relaxed 
since the TDOS is naturally smoothened. Despite sim- 
plicity, this model is well applicable to the experiment [0 . 

In a more complete model, we consider the dissipative 
dynamics of the system (^) within the master equation 
in the Lindblad form Gal, 



so after the null-result measurement, the qubit state, up p = -i[H, p] + V] (■y i /2)([A i p, Aj] + [Ai, pA\]), (5) 



to a phase, becomes ip = A[e^ iEot a o (0)\0) 
where the normalization A = [|ct (0)| 2 + |ai(i)| 2 ]~ 1//2 . 
Generically, A > 1. This corresponds to an increase of 
the component a®, even though the state |0) is fully dis- 
connected. Notice that this result coincides with Eq. (|l|) 
if the term e~ Tt is replaced by |«i | 2 

The integration in Eq. (||) can be formally done as a 
sum over residues e n , the exact eigenvalues of the Hamil- 
tonian (|^). Qualitatively, we can characterize the spec- 
trum of the right well by the average energy spacing A, 
average tunneling amplitude / (the r.m.s. of \f n \ at ener- 
gies near E\), and the total energy bandwidth A ^> A, / 
[see Fig. |l|(c)]. Then at time t < A -1 the contributions of 
different residues add nearly in phase, and |«i| 2 changes well, 
quadratically in t. At t > A -1 the resonant processes of 
return from the right well become important; the form 
of ct\ (t) differs qualitatively depending on the number of 
strongly coupled levels (~ //A) and their exact position. 



where p is the density matrix (DM), Aj are the decoher- 
ence operators, and 7$ are the corresponding decoherence 
rates. We will specifically consider the cases of phase 
noise between the wells, Ao = X]nl n )( n l ( c ^- j0§)' as 
well as the incoherent transitions up and down the lad- 
der of levels in the right well, Ai = A^ = £„ \ n )( n + 1| 
(cf. Ref. 0). 

We have found an exact real-time analytical solution 
of Eq. (§) in the case of uniformly coupled equidis- 
tant states in the right well with infinite bandwidth: 
/„ = /, E n = 71 A, A — * 00. The solution is con- 
structed using the momentum representation in the right 



(2-r l/2 z„ 



Jn<p 



Then the Hamilto- 



nian Hr becomes a differential operator, H R ip(<p) = 
—iA.d c / > 'ijj((f>), the tunneling operator picks <f> = since 
T J ip((p)\(j)) d(j) = |1) v / 27r/-0(O), and the decoherence op- 



erators are diagonal: Aq = 1_r, Ai 
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Then, e.g., the off-diagonal component pio of the qubit 
DM can be found from the equations 

pio = i(E - E x )p w - igb o (0), (6) 
bott) = (iE -'y/2)b (<f>)-Ad^bo(<f>)-ip 10 g8(<f>),('f) 

where 7 = 70 + 71 + 72 is the net dephasing rate, g = 
v2nf, and bo(<f) encodes the off-diagonal components of 
the DM between the level |0) and the right-well levels. 
Eq. (^) is the first-order quasi-linear partial differential 
equation (PDE); it can be integrated in quadratures for 
any form of p±o = pio(t). The solution describes the 
chiral propagation of the decaying amplitude &o around 
the circle from <f> = 0. As a result, before a full turn is 
completed, for t < 2ir/ A, the amplitude does not return 
back to the left well. For this time interval, the only 
effect on the component poi is a relaxation with the rate 
r/2, independent of the dephasing 7: 

Pw (t)=p w (0y( E °- E ^e- rt / 2 , r = 2vr/ 2 /A. (8) 

Similarly, the evolution of the qubit DM component p\\ 
is coupled with the right-well DM components pn(4>, <fi') 
and components &i(</>) involving level |1) and right-well 
levels. Solving in turn the PDEs for pr(4>, 4>') and b\ ((f)), 
we have obtained the self-consistent equation for p\\. 
Again, over the same time interval, there is no return 
tunneling from the right well, and the evolution of p\\ is 
not affected by decoherence described by operators Ao,].^, 
so that pix(t) — pn(G) e~ rt , t < 2ir/A. Evolution of poo 
is trivial: p 00 {t) = p 00 (0). 

With probability poo(t) + pn(t), an ideal projective 
measurement at time t will show that the system has not 
decayed to the right well. In this case the density matrix 
needs to be changed (collapsed) as || |l3f 
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(9) 



It is easy to see that if the system originated in a pure 
state, pl(0) — ipL(0)ip L (0), the density matrix after such 
a measurement at time t < 2it/A also describes a pure 
state, which exactly corresponds to the formula ([!]). We 
emphasize that the absence of return tunneling from the 
right well is both necessary and sufficient for the qubit 
purity in our model. 

The obtained analytical solution of the ideal case can 
be now used as a starting point of the perturbation theory 
for situations more realistic experimentally. In particu- 
lar, a weak non-linearity in the right-well spectrum, e.g., 
E n = A (n + (3n 2 ), f3 <C 1, in the phase representation 
corresponds to a dispersive term 8Hr — —f3Ad 2 /d<f) 2 . 
The analog of Eq. (Q) would then include not only prop- 
agation but also dispersion of the wave packet, and 
the qubit decoherence due to reverse tunneling proce- 
ses may start earlier. The effect is exponentially small 
for 2tt - tA > 5(f), where 5<j) ~ 2{t(3A) 1 ^ 2 is the r.m.s. 
width of the packet. A similar dispersive effect results 
from random level spacing in the right well, or due to 




Yo/A=l., Yi/A=0.8, y 2 /A=1.2 
f/A=0.5, SE/A=-0.2, N=42 



FIG. 2: (Color online) Simulation results for the evolution 
projected onto the left well beginning with the state with 
(a x ) = 1; parameters as shown. Thick lines represent the pu- 
rity P(t), thin lines represent the longitudinal qubit polariza- 
tion — (o z ). Solid lines: equidistant spectrum in the right well, 
E„ — SE + nA; numerical results are almost on top of the an- 
alytical result with TV — > 00 (dotted lines) . The broad minima 
at t > 2n/A and t > 4ir/A correspond to arrival of incoher- 
ent echoes from the right well. Dashed lines: data with added 
spectrum non-linearity, j3 = 0.02 (see text). Dot-dashed lines: 
random level spacing, E n ^.\ — E n — (0.5-1-^^)^, where £n are 
uniformly distributed random numbers, < f n < 1. 

phase noise in the right well. With many levels in the 
right well effectively coupled, we expect these effects to 
be weak at sufficiently early times, as long as the corre- 
sponding phase broadening 5(f) is small, 8<f> <ti 1. 

A different sort of perturbation results if the tunnel- 
ing amplitudes /„ are not equal to each other, or if the 
number of states in the right well is finite. In this case 
the transformed tunneling Hamiltonian would not corre- 
spond to (5-function in phase space but instead couple to a 
finite range of phases, |(/>| < £0 ~ A/A. As a result, some 
reverse tunneling from the right well back to the qubit 
may start early. However, the associated decoherence is 
not expected to be significant as long as 8(p <C 27r — tA. 
Additionally, the decay in the model (Jsj) would be expo- 
nential only at t ^> A -1 , which leads to an additional 
prefactor as in Eq. (||). We also note that the inelastic 
escape to levels decoupled from |1) also reduces the re- 
turn tunneling probability and extends the time interval 
of qubit coherence. 

We illustrate these arguments by a numerical simula- 
tions of Eq. (||) in Fig. ||: thick lines represent the purity 
P(t) = [(fiu -Poo) 2 +4|poi| 2 ] 1/ 7(A)o +pn), while thin 
lines represent the diagonal component of qubit polariza- 
tion, -(cr z ) = (poo-pn)/(poo + pn)- Even with not very 
large parameter A/ A = 40 (which corresponds to total 
of N — 42 energy levels), the results of numerical sim- 
ulation agree almost perfectly with the analytical result 
shown by dotted lines. In agreement with our arguments, 
neither weak spectrum non-linearity (dashed lines) nor 
randomized level spacing reduce the qubit coherence for 
sufficiently early evolution time. 

Let us briefly discuss what happens with Eq. ([!]) when 
the state |0) can also tunnel into the right well with the 
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rate Tq. The tunneling would always produce a matrix 
element hio between the levels |1) and |0). In the sim- 
plest case of no return processes, such a matrix element 
would be purely imaginary (anti-Hcrmitian) , with the 
magnitude \h w \ < (r r) 1 / 2 /2 @, where the equality 
is achieved when the tunneling matrix elements from the 
two states are all proportional to each other, To n = (Tin- 
The corresponding dynamics is not trivial, but the level 
mixing is small if \h\o\ -C \E\ — Eq\. Even in this case 
there may be a significant effect on level decay rate. How- 
ever, if r + Ti <C \Ei — E \, then Eq. (P can be simply 
replaced by the result expected from the quantum Bayes 
rule)?], Q: the term a (0) should be substituted with 
a (0)exp(-rV/2). 

So far we have assumed an ideal orthodox measure- 
ment at time t, which distinguishes between the left and 
right wells, but does not affect the states in the left 
well. For a superconducting phase qubit such a measure- 
ment can be technically realized by biasing the measure- 
ment SQUID at time t to the point at which the SQUID 
switches to the finite-voltage state only for a right-well 
fluxflo||. Due to strong nonlinearity of such a detector, 
in the case of no switching the back-action onto the left- 
well qubit states can be made practically negligible. It is 
important to mention that a linear detector would neces- 
sarily disturb the qubit states because for a phase qubit 
the "distance" between the wells is comparable to the 
"width" of the qubit well. 

Strictly speaking, the results discussed in this work as- 
sume measurement only once at time t. However, there 
is a sense in which our results describe qubit evolution 
before the decay. Our qubit "ages" in the process of no 
decay [as in Eq. (P], unlike the case of a radioactive 
atom which remains "as new" before the decay actually 
happens. For such an interpretation we necessarily need 
to consider repeated (or continuous) measurements with 
time resolution better than T -1 , and it is important that 
presence or absence of extra measurements within the 
time interval t should not affect the non-decayed qubit 
state at time t. Clearly, this is not the case if ortho- 
dox measurements are repeated too frequently with time 
interval At shorter than the scale of the quantum Zeno 
effect, At < A -1 . However, in the regime of an expo- 
nential decay with At ^> A -1 [e.g., as in Eq. (|J) with 
Z « 1] an extra measurement has no effect: a compo- 
sition of two evolutions (P with durations t\ and t% is 
the same as a similar evolution with duration t\ + 1%. In 
actual experiment the measurement SQUID can monitor 
the decay continuously, and then At corresponds to the 
intrinsic time resolution of the detector. In this case the 
frequent partial collapses can be replaced by introduction 
of the interwell phase noise [Eq. (JSJ)] with 70 = 1/At. 
Our results indicate that this does not lead to significant 
qubit dephasing even in the case of good time resolution 
At <C r _1 , as long as the conditions of exponential de- 
cay are well satisfied. We conclude that in the regime of 
exponential decay, Eq. dll) can really be interpreted as 



actual qubit evolution in time before decay. 

In conclusion, we have analyzed quantum dynamics of 
a model with two-well structure resembling the experi- 
ment ||, with the (qubit) states in the left well nearly 
coherent, while those in the right well and the transition 
Hamiltonian are subject to decoherence. The analytical 
solution of the master equation (|^) obtained for infinitely 
wide right-well spectrum with equal level spacing A gives 
pure qubit subspace for t < 2ir / A, independent of the de- 
coherence rates 70,1,2, see Eq. (||). This property of co- 
herence preservation over finite time interval remains in 
effect in a perturbed system where the solvability condi- 
tions are only approximate. We have identified a regime 
where the quantum evolution during tunneling can be ex- 
perimentally accessed via a repeated "stroboscopic" mea- 
surement or a continuous "weak" measurement. In this 
regime the qubit state will remain pure in spite of the 
phase noise associated with the measurement. 
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